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We investigate high density state of SU(2) QCD in the case of the Wilson fermions with Iwasaki improved 
action. Based on the Woodbery formula, the ratio of fermion determinants is evaluated at each step of the 
Metropolis link update. At P — 0.7, for k = 0.150, and « = 0.175, we calculate thermodynamical quantities, such 
as, baryon number density, Polyakov line, and the energy density of gluon sector with chemical potential fj,=0 to 
0.9 on the 4 3 x 8 lattice. The susceptibility of the Polyakov line to the chemical potential shows peak at around 
fj, — 0.7 which indicates the deconfinement phase transition. Behavior of the meson and diquark propagators with 
finite chemical potential are also investigated at both side of the peak. . 
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1. Introduction 

High density state of the strongly interacting 
matter is a very attracting problem in the con- 
nection to the recent experiments of the Rela- 
tivistic Heavy Ion Collider at BNL jlj]. Physics 
of extremely high density state is also important 
for the understanding the astrophysical objects 
such as neutron star. Usually in the statistical 
physics, high density state is investigated with 
use of the ground canonical ensemble with chem- 
ical potential which is introduced as a Lagrange 
Multiplier for the conserved charge. However, as 
is well known, chemical potential makes action 
of the QCD complex and causes a problem in the 
numerical convergence of simulation based on the 
monte calro method 

SU(2) gauge theory is the simplest non-abelian 
gauge theory and holds the pseud-real property 
which keeps action real even with chemical po- 
tential. Though the real QCD is not SU(2) but 
SU(3), SU(2) gauge theory is expected to play an 
role as a QCD-like theory which informs us im- 
portant features of the non-abelian field theory 
j3|. SU(2) gauge theory as a residual interaction 
in the Color Super Conducting state is also one 
of the promising possibility and the high density 
state of SU(2) gauge theory can be realized B. 



We calculate thermodynamical quantities as 
functions of chemical potential. Ground canon- 
ical ensemble states are generated through up- 
date based on the Woodbery formula Q . Iwasaki 
improved action is adopted for the gauge action. 
Though the improved action seems to stabilize 
the simulation better than the naive Wilson ac- 
tion Q, our simulation still shows instability at 
around = 1.0. However, the susceptibility of 
the Polyakov line shows peak at around 

/i = 0.7 which suggests the deconfinement phase 
transition. Therefore, our simulation can corre- 
spond to the state around across the phase tran- 
sition. 

2. Thermodynamical quantities 

With use of Wilson fermions and Iwasaki 
improved action we calculate thermodynamical 
quantities, such as Polyakov line, gluon energy 
density and baryon number density with 4 3 x 8 
lattice. The expectation value of the baryon num- 
ber density is given by, 

<n> =w s i logZ (1) 

where V s is the spatial volume N x x N y x N z . 
Energy density, e, is given as, 
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Figure 1. Polyakov line, gluon energy density, and 
baryon number density as a function of chemical 
potential /i. 
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Figure 2. Susceptibility of the Polyakov line as 
a function of chemical Potential fj,. The circle 
plots and diamond plots correspond to the k = 
0.150 and k = 0.175, respectively. The peak point 
seems to depend on the hopping parameter k 



The derivative of the partition function is com- 
posed of two parts, 

(log Z)' = - J VUV^(-S' G -S' F )e^- SG - s -\(3) 

where Sq and Sf are gluon action and fermion 
action, respectively. We denote the contribution 
of the gluon action part (first term in the r.h.s. 
of eq.(3)) by gluon energy density. Figure 1 dis- 
plays Polyakov line, gluon energy density and 
baryon number density as functions of /i where 
(3 = 0.7 and n = 0.150, respectively. All quanti- 
ties start to have non-zero value at about /i = 0.4 
and rise up with chemical potential in the region 
0.4 < /i < 0.8. None of them is order parame- 
ter of the phase transition in the exact sense, and 
since the lattice size is small, no sharp change is 
seen. However, growing up of the these quanti- 
ties indicates that quarks and gluons become free 
from the confinement force at finite chemical po- 
tential about 0.4 < fi < 0.8. 

Figure 2 shows susceptibility of the Polyakov 
line as the functions of the chemical potential 
fi. The peak appears at around /i = 0.7 which 



suggests that our simulation crosses the decon- 
fincmcnt phase transition. However, peak point 
seems to depend on the hopping parameter k. 
Larger hopping parameter seems to give lower 
phase transition /ic- 

3. Propagators 

Figure 3 shows our preliminary results of the 
propagator of mesons and diquarks with k = 
0.175. In the lower figure, at fi — 0.0, pion 
and scalar diquark(b5b), and aO and pseudscalar 
diquark(blb) are degenerate each other, as ex- 
pected. However, such degeneration of meson and 
diquark disappears with finite chemical potential. 
Because diquark possesses definite baryon num- 
ber, the effect of the chemical potential on di- 
quark and anti-diquark are in the opposite sign. 
Therefore, the symmetry of the propagator in 
time direction and anti-time direction is broken 
and asymmetry in nt appears(blb and b5b in the 
lower figure of Fig. 3). On the other hand, mesons, 
which have vanishing baryon number, keep sym- 
metric property in nt direction 0. According to 
figure 2, /i = 0.8 with k = 0.175 is in the region 
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Figure 3. Preliminary results of the propagator 
of the meson and diquark with k = 0.175. The 
upper figure and lower figure stands for the fi = 
0.8 and ^ = 0.0, respectively. If we look upon 
the peak in figure 2 as the deconfinement phase 
transition point, those figures correspond to the 
mesons diquarks in the state over and under the 
phase transition 



upper side than the peak point. 

4. Concluding Remarks 

SU(2) finite density state is investigated with 
Wilson fermion and Iwasaki improved action. 
Peak of the susceptibility appears at around fj, = 
0.7, but it seems to depend on the k. It means 
physical size of the lattice spacing will change 
with k and careful treatment will be required in 
the procedure of the chiral limit. The propaga- 
tors of the mesons and diquark state (Baryon in 
SU(2)) have also discussed JtJ. Detailed analyses 



for both meson and diquark will appear some- 
where else. 
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